Periodically-dressed Bose-Einstein condensates: a superfluid with an anisotropic and 

variable critical velocity 
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Two intersecting laser beams can produce a spatially-periodic coupling between two components 
of an atomic gas and thereby modify the dispersion relation of the gas according to a dressed- 
state formalism. Properties of a Bose-Einstein condensate of such a gas are strongly affected by this 
modification. A Bogoliubov transformation is presented which accounts for interparticle interactions 
to obtain the quasiparticle excitation spectrum in such a condensate. The Landau critical velocity is 
found to be anisotropic and can be widely tuned by varying properties of the dressing laser beams. 



In numerous physical systems, particles which are con- 
fined to a medium can be treated as free particles whose 
properties are modified by the medium. For example, 
the dispersion relation of electrons in a periodic solid is 
changed from the free-particle relation to a band struc- 
ture which is periodic in quasi-momentum space. Such 
electrons are then treated as particles with properties 
which can be purposely engineered by modifying the sur- 
rounding periodic structure. In this Letter, we consider 
similarly the task of engineering novel behaviour in a 
quantum-degenerate atomic gas by placing the gas in a 
properly constructed periodic medium. 

A periodic potential for an atomic gas can be produced 
by intersecting two or more laser beams. A polarizable 
atomic gas illuminated by two intersecting off-resonant 
laser beams with identical polarization and frequency but 
different wavevectors experiences a spatially-periodic po- 
tential proportional to the light intensity, i.e. the atoms 
reside in a crystalline medium made of light. Atoms 
in such media have been studied in the non-degenerate 
and quantum degenerate regimes. A small fre- 
quency difference to = LO1—UJ2 can be introduced between 
the laser beams to induce resonant Bragg transitions be- 
tween identical internal but different momentum states. 
Such Bragg transitions have been used to probe Bose- 
Einstein condensates . 

In this Letter, we present a different scheme for en- 
gineering properties of an ultra-cold gas. Rather than 
considering a spatially periodic potential, we consider a 
spatially-periodic coupling between two internal states of 
an atomic gas, i.e. we consider a periodic medium con- 
tructed by laser beams which effect Raman rather than 
Bragg transitions. An atomic gas in this medium is char- 
acterized by an easily tunable dressed-state dispersion 
relation. In creating a Bose-Einstein condensate of such 
a periodically-dressed atomic gas, one may explore how 
the nature of free single-particle excitations can modify 
macroscopic properties of a quantum fluid. In particular, 
we develop a Bogoliubov transformation which defines 
the quasi-particle excitation spectrum, and show that 
this novel quantum fluid has a tunable and anisotropic 
superfluid critical velocity. 

Let us consider a uniform dilute gas composed of atoms 



of mass m with two internal ground states, |o) and \b), 
separated by an energy htoo, and an excited state |e) (see 
Fig. . This gas is exposed to two laser beams (labeled 
1 and 2), with wavevectors ki and k2 and frequencies wi 
and LU2, respectively. The beams are polarized so that 
beam 1 connects states \a) and |e) and beam 2 connects 
states 1 6) and |e), while both beams have a large de- 
tuning A from these allowed dipole transitions. These 
laser beams can induce Raman transitions between the 
two internal ground states. A transition from state a) 
to b) : results in a momentum kick of Kk = ?i(ki — k2). 
Such coupling can be represented in position space as 
a spatially-periodic coupling between the two internal 
states, proportional to e''' ''|6) (a| + e~*'^ '"|a) 
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FIG. 1. Engineering propeties of a periodically-dressed 
atomic gas. (a) Laser beams of frequency lji and 0)2 may 
induce Raman transitions between internal states |a) and \b). 
S = (iui —LJ2)~u}o is the detuning from the Raman resonance, 
(b) Such a Raman transition imparts a momentum transfer 
of ?ik = ?!.(ki — k2), where ki and k2 are the wavevectors of 
the Raman coupling lasers. 

Using second quantized notation, a gas which is con- 
tinuously illuminated by these beams can be described 
by the dressed-state Hamiltonian 



a' a„ + o'b„ 

2m 2/'i^ V2m 2''i^ 



E («4^q+k/2Ciaq-k/2 + e*cla]^_i^/2C26q+k/2) (1) 
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Here Oq and bq (ajj and b^) are the annihilation (cre- 
ation) operators for atoms with wavevector q in inter- 
nal states I a) or respectively. The operators ci and 
C2 {c\ and cj), which are photon annihilation (creation) 
operators for photons in beams 1 and 2, can be re- 
placed with c-numbers y^Ni and V N2, respectively, in 
the limit that the photon number is large. We then define 
n = 2e^/{Ni + 1){N2 + l)/h as the (real) two-photon 
Rabi frequency. 

Now define single-atom states |aq) = |A^i -1-1; A^2; a, 
k/2) and |6q) = |^i;iV2 + l;&,q+ k/2) which are con- 
nected by a Raman transition. The notation indicates 
that for state |aq) there are A^i-f 1 photons in beam 1 and 
N2 photons in beam 2, and the atom is in state |a) with 
momentum tiq + k/2, and similar for state |&q). These 
two states have the same total momentum Q , which we 
define as hq by a proper choice of reference frame. The 
unperturbed energies (e — s- 0) of these states are 



(2) 

where S = {uji — 0^2) — ti^o is the detuning from the Raman 
resonance. 

Diagonalizing the Hamiltonian in the subspace of 
states |aq) and |5q) yields the dressed-state atomic eigen- 
states |±q) with energies 
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huj„ = I q~ 



(4) 



and creation operators TTq (for state 
state |— q)) defined as 



"q 

Mq 



cos6'q/2 smdq/2 
-sin0q/2 cos(?q/2 



-q)) and ^q (for 



(5) 



The dressed states are derived from the bare states |aq) 
and |6q) by a rotation in the two-dimensional subspace 
by the angle —0q/2 where tanSq = il/{d + 'hq- k/m). 

The dressed-state dispersion relation (Eq. ^ is shown 
in Fig. ||. The motional ground state occurs in the lower 
dressed state at a wavevector Q, which is near either 
±k/2 depending on the sign of S. 
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FIG. 2. Dressed-state dispersion relation: single parti- 
cle energies hw^ vs. (normalized) wavevector q/k. The 
excitation momentum ?iq and the momentum transfer 
Kk = /i(ki — k2) are chosen to be colinear. Energies are 
scaled by = Ti^ /2m,. For the traces shown, 

Ml/Ef^l2 = 2, and h5/Ek/2 ~ 3 (solid lines), h5/Ei./2 = 
(dashed lines), or Ti5/ Ef^/2 ~ —1 (dotted lines). For small \S\ 
(e.g. h5/Ek/2 ~ —1), the dispersion relation acquires a sec- 
ondary minimum. As S changes sign, the momentum of the 
lowest energy state changes discontinuously. For 5 = 0, the 
motional ground state is degenerate. 

Let us now consider the effects of Raman coupling on a 
two-component Bose-Einstein condensate, which is now 
formed of a macroscopic population of atoms in the low- 
est energy state |— q), i.e. the condensate is formed of 
atoms in a coherent superposition of atomic states with 
differing momenta laQ-k/2) ^-nd |6Q+k/2)- This conden- 
sate has a two-branch excitation spectrum, reflecting the 
presence of two internal states (and, correspondingly, two 
distinct dressed-state levels). As with a scalar Bose- 
Einstein condensate, weak repulsive interactions should 
yield phonon-like excitations at low energies and a free- 
particle-like dispersion relation at high energies. How- 
ever, unlike for a scalar condensate, we expect properties 
of a periodically-dressed Bose-Einstein condensate to re- 
flect the anisotropy of the dressed-state dispersion rela- 
tion. 

Let us now treat explicitly the effects of atomic colli- 
sions by writing the many-body Hamiltonian as 



H = ^ (nWq /i]jpiq -I- ?iW+7rt TTq) + T^j^t 

q 



(6) 



Considering only elastic binary collisions (which conserve 
the number of atoms in each of the internal states) char- 
acterized by identical s-wave scattering lengths a, we may 
write the interaction Hamiltonian Ti.[^i as 



^int-fE(^ 



N) 



(7) 



Here g = (47r ?i^a/m) x V^^, V is the volume occupied 
by the gas, N is the number of atoms in the gas, and riq, 
the spatial Fourier transform of the density operator, is 
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Note that, since the rotation matrix TZ (defined in ac- 
cordance with Eq. 1^) is not generally diagonal, elastic 
collisions do not generally conserve the number of atoms 
in the upper and the lower dressed states, respectively. 

We now make use of the Bogoliubov approximation 
1^ in which we assume a macroscopic population of A^o 
atoms in the lowest energy state of wavevector Q, and 
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substitute /iq = /iq — \ 
the 4-coniponent vectors 



To proceed, we consider potential /i = /ix 2.5 kHz corresponds to a condensate 
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(10) 



The Bose commutation relations of the dressed-state an- 
nihilation and creation operators can be expressed as 



Isolating terms of order N and iV, we 



may now approximate the Hamiltonian as 



, 2 



(11) 



where H is a 4 x 
The diagonal matrix £ij 



4 matrix of the form Hij = £ij 



q where £q 



has entries £q , — fZq, , and 
The chemical potential 



-OJ, 



Q, 



is given hy n = gN , and we define 



cos Aq 
-i cos A_c 

sin Aq 
-i sin A_o 



(12) 



where Aq = (6'Q+q - 6»q)/2. 

The quasi-particle energies and their creation and 
annihilation operators are found by diagonalizing the 
matrix H. Given the invertible matrix M for which 
MHM^^ = H is diagonal, allowing v to transform as 
a column vector (vi = MijVj) and w as a row vector 
{wj = WiM~j^), we find the Bose commutation relations 
[vi^Wj] — Sij to be preserved for the quasi-particle cre- 
ation and annihilation operators (/iq, /i^, 7fq, tt^) defined 
by the elements of v and w. One finds the diagonal ele- 



ments of ?i iJ to be Cjr 



and 



^Q+q' ~"^Q-q' '^Q+q' '^^^^ ^"^Q-q 

which define the lower and upper quasi-particle excita- 
tion energies at wavevector Q ± q. 

Figure ^ shows the quasiparticle spectrum calculated 
for excitations parallel to the Raman transition momen- 
tum transfer k. In choosing parameters for this calcula- 
tion, we have in mind an experimentally convenient re- 
alization with a Bose-Einstein condensate of ^^Rb. One 
may choose the internal ground hyperfine states \a) = 
\F = l,mF = —1) and \b) = \F — 2,mp — 1) which can 
be connected by a two-photon Raman transition. Choos- 
ing these states has the benefit that the Raman transition 
frequency is insensitive to magnetic field fluctuations, 
both hyperfine states are magnetically trappable, and, 
specifically in ®^Rb, inelastic collisions are scarce. Fur- 
thermore, the scattering lengths for all elastic collisions 
are nearly the same, justifying the assumption made in 
Eq. 1^. We consider the case of counter-propagating Ra- 
man laser beams which are nearly resonant with the D2 
transition (fc ~ 47r/A where A — 780 nm). The chemical 



density of 3 x 10 cm 
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FIG. 3. Bogoliubov quasi-particle dispersion relation of 
a periodically-dressed Bose-Einstein condensate. We con- 
sider a gas of *^Rb exposed to counter-propagating Raman 
beams near the *^Rb principal transition, with hS = ~Ek/2, 
Ml = 2Ek/2 and fi = 0.61i5fe/2 = hx 2.5 kHz. Excitations 
parallel to k are considered. Lower and upper quasi-particle 
excitation branches are shown (solid lines) as a function of 
the wavevector q of the excitation. The quasi-particle en- 
ergy spectrum is higher than the free dressed-state dispersion 
relation (dashed lines). Dotted lines show the quasi-particle 
dispersion relation in the case of no Raman coupling (n — > 0). 
The effect of Raman coupling is to induce a level anti-crossing. 

As shown in the figure, the quasi-particle energies are 
higher than the free dressed-state energies due to the 
repulsive interactions between atoms. Comparing the 
quasi-particle spectrum to that for the two-component 
condensate in the absence of Raman coupling (il 0), 
one sees that the effect of the dressing lasers is to intro- 
duce avoided level-crossings to the spectrum. The spec- 
trum for quasi-particle excitations near the condensate 
momentum Q (i.e. wavevectors Q -f q for small q) is lin- 
ear, describing phonon-like excitations. For excitations 
parallel to the direction of momentum transfer and mak- 
ing small q approximations to the Hamiltonian, we find 
the lower excitation spectrum t o have the limiting value 
^'^Q+q ~ c*hq where c* = a/ fj,/m* is the Bogoliubov 
speed of sound corresponding to an effective mass m* 
determined by the curvature of the dressed-state disper- 
sion relation at its minimum. We note further that the 
complex matrix H which appears in the Hamiltonian has 
two positive and two negative real eigenvalues |l^ , as is 
required for the stability of the condensate. 

Finally, we consider the implications of the dressed- 
state dispersion relation for the superfluidity of the 
periodically-dressed Bose-Einstein condensate. An ex- 
planation for the dissipation-less flow of a superfluid be- 
low a critical velocity was provided by Landau, who 
used kinematic arguments to define — imiiE{q)/hq 
where E(q) is the quasi-particle excitation energy at 
wavevector q For weakly interacting scalar Bose- 

Einstein condensates, the Bogoliubov quasi-particle dis- 
persion relation gives a Landau critical velocity = c 
which is equal to the speed of sound c = -y/ n/m, and 
the onset of scattering of microscopic impurities thus oc- 
curs by the scattering of phonons. This contrasts with 
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the behaviour of superfluid ^He in which excitations are 
phonon-Uke at low momenta, while at higher momenta 
there exists a secondary minimum in the excitation spec- 
trum corresponding to rotons fl^ j. The Landau critical 
velocity in that case is set by the roton minimum [ p^ . 

The presence of a secondary minimum in the disper- 
sion relation of periodically-dressed Bosc-Einstcin con- 
densates suggests an analogy to superfluid ^He and a 
reduction of the superfluid velocity below the speed of 
sound. Figure ^ shows the Landau critical velocity cal- 
culated for impurity velocities parallel to the Raman mo- 
mentum transfer k, as the Raman laser detuning S is var- 
ied. For large values of \S\, the critical velocity is equal 
to the speed of sound \Jiilm as for a single-component 
Bose-Einstein condensate. As is lowered, the critical 
velocity becomes anisotropic. In one direction, is dra- 
matically lowered due to scattering at the secondary min- 
imum (an "artificial roton") in the dispersion-relation. 
We may approximate the "artificial roton" minimum as 
occurring at energy 'h\i\ and quasi-particle momentum 
Kk, obtaining \vi\ — |<5/fc| in this regime. Thus, the su- 
perfluid properties can be controlled by varying the de- 
tuning b or by varying k by changing the relative angle 
between the Raman laser beams. In the other direction, 
— c* — ^ iijra* as determined by phonon scattering. 
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FIG. 4. Landau critical velocity in a periodically-dressed 
Bose-Einstein condensate. Velocities aligned with (positive 
vl) or counter to (negative vl) the Raman momentum trans- 
fer k are considered, and we take fiQ. = 2Ei^/2- At large 
detunings \5\, vl has the same magnitude for flow in both 
directions, with a valu e approaching the Bogoliubov speed 
of sound c — -y/ jj./m — 3.3 mm/s for a ^'^Rb condensate 
at the density of 3 x 10^'' cm^'^. At smaller detunings, an 
anisotropy in vl develops. For the regimes {S > 0, vl > 0} 
and {5 < 0, Vl < 0}, Vl is determined by the speed of Bogoli- 
ubov sound modified by the effective mass m* > m. In the 
regimes {S > 0,vl < 0} and {5 < 0,vl > 0}, the magnitude of 
Vl is lowered due to the secondary minimum in the dispersion 
relation. 

One may also vary the intensity of the Raman beams, 
thereby changing f2. An important impact of changing 
the Raman coupling strength is not only altering the Lan- 
dau critical velocity (changing both the effective mass 
TO* and the position of the secondary minimum), but 
also controlling the degree to which quasi-particles in 
the lower or upper excitation branches can be created by 
scattering off an obstacle. That is, the Landau criterion 



determines the onset of dissipation for a moving super- 
fluid, but does not describe the strength of such dissipa- 
tion. A large Raman coupling would enhance scattering 
into secondary minimum of the lower excitation branch, 
while in the limit — > 0, this scattering rate clearly 
vanishes. A detailed calculation of this dissipation rate 
(similar to Refs. |I|-|l|]) will be given elsewhere. 

In summary, we have described a means of engineering 
a novel quantum fluid composed of dressed-state atoms 
in a spatially-periodic Raman coupling medium. This 
quantum fluid should be amenable to study using cur- 
rent methods for probing ultracold atomic gases. The 
quasi-particle dispersion relation can be probed by study- 
ing collective excitations (i.e. density and magnetization 
modulations) at various length scales (see reviews in 
JlTf). Aspects of superfluidity in gaseous Bose-Einstein 
condensates have been investigated in recent experiments 
||l8|,^, and similar methods can be used to investigate 
the superfluidity of a periodically-dressed condensate. 
For example, the anisotropic Landau critical velocity can 
be probed by studying the scattering of an impurity gas 
which propagates through the fluid ||l5| . Further the- 
oretical work should address a number of issues. For 
instance, the current treatment describes the quantum 
fluid exclusively in momentum space. Experimentally, 
such a fluid would be produced within a trapped volume. 
Excitations with wavelengths much smaller than the size 
of the fluid (as probed with Bragg scattering |^) should 
be well described by our treatment, while the descrip- 
tion of longer wavelength excitations will require further 
theoretical considerations. Another interesting consid- 
eration is the behaviour of the quantum fluid with the 
Raman lasers tuned precisely to the Raman resonance, 
5 = 0, at which point the motional ground state is dou- 
bly degenerate. We suspect that in the presence of an 
external trapping potential, such a quantum fluid would 
be described as evolving in a double- well potential in mo- 
mentum space. Finally, our scheme may also be extended 
to condensates with more than two internal states, such 
as spinor Bose-Einstein condensates |l9| , giving greater 
flexibility in engineering properties of gaseous quantum 
fluids. 
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